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ABSTRACT 

We analyze X-ray spectra of 43 Palomar-Green quasars observed with XMM-Newton in order to 
investigate their mean Fe K line profile and its dependence on physical properties. The continuum 
spectra of 39 objects are well reproduced by a model consisting of a power law and a blackbody 
modified by Galactic absorption. The spectra of the remaining four objects require an additional 
power-law component absorbed with a column density of ~ 10 cm . A feature resembling an 
emission line at 6.4 keV, identified with an Fe K line, is detected in 33 objects. Approximately half 
of the sample show an absorption feature around 0.65-0.95 keV, which is due to absorption lines 
and edges of O VII and O VIII. We fit the entire sample simultaneously to derive average Fe line 
parameters by assuming a common Fe line shape. The Fe line is relatively narrow (a = 0.36 keV), 
with a center energy of 6.48 keV and a mean equivalent width (EW) of 248 eV. By combining black 
hole masses estimated from the virial method and bolometric luminosities derived from full spectral 
energy distributions, we examine the dependence of the Fe K line profile on Eddington ratio. As the 
Eddington ratio increases, the line becomes systematically stronger (EW = 130 to 280 eV), broader 
(a = 0.1 to 0.7keV), and peaks at higher energies (6.4 to 6.8keV). This result suggests that the 
accretion rate onto the black hole directly influences the geometrical structure and ionization state of 
the accretion disk. We also examine a two-component model consisting of a Gaussian and a diskline 
to constrain the intensity of the broad line. The mean equivalent widths are « 70 — 180 eV for the four 
Eddington ratio groups, although the standard deviations in each group are very large. This suggests 
that the broad line is not ubiquitous. 

Subject headings: accretion, accretion disks — galaxies: active — galaxies: nuclei — galaxies: Seyfert 
— quasars: general — X-rays: general 



1. INTRODUCTION 

Fluorescent Fe K lines are common in the X-ray spectra 
of accreting black holes, ranging from X-ray binaries to 
active galactic nuclei (AGNs), and provide one of the 
best probes for studying accretion disks. If the Fe K 
line is emitted from the inner part of the accretion disk, 
it becomes broad and asymmetric due to both Doppler 
shift and gravitational redshift (Fabian et al. 1989). The 
best example of a broad Fe K line is seen in the Seyfert 
1 galaxy MCG-6-30-15 (Tanaka et al. 1995). Broad Fe 
K lines have been observed in many other sources and 
are thought to be very common in the X-ray spectra of 
Seyfert 1 nuclei (Nandra et al. 1997a). 

Recent observations of AGNs with XMM-Newton show 
that the broad Fe K line is not as common as previously 
believed. While a few Seyfert 1 galaxies (e.g., MCG-6- 
30-15, Mrk 205, Mrk 509) indeed show an unambigu- 
ous broad line, Fe lines in other galaxies are dominated 
by a relatively narrow feature. From an analysis of 53 
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type 1 AGNs observed with XMM-Newton, Page et al. 
(2004a) showed that a broad Fe K line (a > 0.1 keV) is 
seen in only 13 sources, while the remaining 40 sources 
have a narrow line (a < 0.1 keV) or no lines. Some of 
the 13 sources show a narrow line as well as the broad 
line. In a similar study, Yaqoob & Padmanabhan (2004) 
reported the results of 18 observations of 15 Seyfert 1 
galaxies observed with the Chandra High Energy Grat- 
ing. They measured the width of the line core and ob- 
tained a weighted mean of FWHM = 2380± 760 kms -1 , 
which is slightly larger than the instrument resolution 
(FWHM w 1860 kms -1 ). Evidence of underlying broad- 
line emission was also found in at least four sources. 
Based on these recent studies, it is still controversial 
whether or not relativistically broadened Fe K emission 
is truly common in nearby AGNs. 

Streblyanska et al. (2005) derived an average rest- 
frame spectrum of AGNs detected in a 770 ksec XMM- 
Newton observation of the Lockman Hole field. They 
used a sample of 104 X-ray sources with optical redshifts 
measured by Lehmann et al. (2001), analyzing separately 
the type 1 and type 2 subsamples defined by optical spec- 
tra (Schmidt et al. 1998; Lehman et al. 2001). From 
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composite rest-frame spectra generated for the type 1 
and type 2 sources, they found evidence for a broad line 
peaking at a rest-frame energy ^6.4 keV with an equiv- 
alent width (EW) ~560eV and ^460 eV, respectively. 
However, it should be noted that there are some system- 
atic uncertainties in modelling the continuum, since it 
is virtually impossible to distinguish a very broad line 
clearly from the continuum for sources as faint as those 
analyzed by Streblyanska et al. (2005). In addition, 
making composite spectra of faint sources may introduce 
spurious spectral features (Yaqoob 2005). 

Brusa et al. (2005) and Comastri et al. (2007) studied 
average spectra of AGNs detected in the Chandra Deep 
Field North and South. They stacked the X-ray counts 
in the observed frame from spectroscopically identified 
AGNs, using a large number of source spectra within 
sufficiently narrow rcdshift ranges such that the energy 
(rcdshift) spread is negligible. A broad Fc line is seen in 
some of the stacked spectra and is fitted with a relativis- 
ts disk line model peaking at a rest-energy of 6.4 keV. 
Since the fluxes of the objects in the sample are low, this 
study may suffer from the same uncertainties affecting 
the Streblyanska et al. analysis of the Lockman Hole. 

What parameters of the AGN affect the strength or 
profile of the Fe K line? An inverse correlation between 
the equivalent width of Fe K line and X-ray luminosity, 
known as the X-ray Baldwin effect, was first pointed out 
by Iwasawa & Taniguchi (1993). This correlation has 
been confirmed with ASCA (Nandra et al. 1997b) and 
XMM-Newton data (Page et al. 2004a; see also Jimenez- 
Bailon et al. 2005). The dependence of the iron line 
profile with luminosity was investigated by Nandra et 
al. (1997b) using 18 Seyfert Is and 21 quasars observed 
with ASCA. They divided the sample into five luminos- 
ity bins and examined average Fe K line profiles. The 
line shows a very similar profile, which is composed of 
a narrow core and a broad red wing for groups with a 
luminosity below Lx ~ 10 44 ergs s . The intensity of 
the narrow core, however, becomes weak above this lu- 
minosity. The red wing becomes weak and the peak en- 
ergy shifts to higher energy for the luminosity range of 
10 45 < Lx < 10 46 ergs s _1 ; no evidence for line emission 
is observed above Lx ~ 10 46 ergs s _1 . Jimenez-Bailon et 
al. (2005), analyzing a sample of 39 Palomar-Green (PG; 
Schmidt & Green 1983) quasars observed with XMM- 
Newton, found a similar result as that reported by Nan- 
dra et al. (1997b). The luminosity dependence of the 
Fc K line may be related to the ionization level of the 
accretion disk. 

This paper examines the behavior of the Fe K line 
as a function of accretion rate, as parameterized by 
the Eddington ratio, taking advantage of the availabil- 
ity of black hole mass and bolomctric luminosity esti- 
mates. The Eddington ratio is a fundamental parameter 
that may control both the ionization state and geomet- 
ric structure of the accretion disk. If so, we expect the 
profile of the Fe K line, in particular its width and cen- 
tral energy, to vary systematically with Eddington ratio. 
The goal of this study is to investigate this issue. 

This paper is organized as follows. Section 2 describes 
the observations, and Section 3 shows the results of the 
spectral analysis of average spectra and its Eddington 
ratio dependence. We discuss the physical origin of the 
Eddington ratio dependence in Section 4. Our findings 



are summarized in Section 5. 

2. XMM-NEWTON OBSERVATIONS 

2.1. The Sample 

We selected our sample from the 87 low-redshift (z < 
0.5) Palomar-Green quasars in Boroson & Green (1992). 
X-ray data for 45 objects are available in the XMM- 
Newton Science Archive (XSA). Only EPIC-PN data 
were used in the following analysis since the effective 
area of PN is much larger than that of EPIC-MOS. PG 
1001+045 was excluded from the present study because 
no EPIC-PN data are available. PG 1226+023 (3C 273), 
a very bright quasar, was also excluded since it over- 
whelms the X-ray flux in the average spectrum exam- 
ined below. Basic data (name, redshift, central black 
hole mass, and Eddington ratio) for the final sample of 
43 objects are given in Table 1. XMM-Newton results of 
34 objects have been published, and the references are 
also given in Table 1. A log of the observations is given 
in Table 2. 

The black hole masses and Eddington ratios are taken 
from the study of L. C. Ho (2007, in preparation). Ho 
calculates black hole masses employing the virial method 
of Kaspi et al. (2000), using broad H/3 line widths from 
Boroson & Green (1992) and 5100 A continuum lumi- 
nosities from Ncugebauer et al. (1987). His bolomet- 
ric luminosities come from direct integration of the full, 
broad-band (radio to X-rays) spectral energy distribu- 
tion rather than from a single band assuming a bolomet- 
ric correction. All distance-dependent quantities were 
derived assuming the following cosmological parameters: 
O m = 0.3,17a = 0.7, and H = 72 kms Mpc -1 . Edding- 
ton luminosities, computed assuming a composition of 
pure hydrogen, are given by 

4TrGM BH m p c n38 ( M BH \ _i 

^Edd = — = 1-26 x lO' 30 — — ergss \ 

<7 T V M & J 

(1) 

where Mbh is the mass of the black hole, m p is the proton 
mass, and <tt is the Thompson scattering cross-section. 

Figure 1 shows the distributions of redshift, black hole 
mass, and Eddington ratios. The mean redshift is z k 
0.2, but the sample covers a wide range of black hole 
masses (10 6 < M B h/M© < 10 9 ) and Eddington ratios 
(0.05 < iboi/^Edd < 4.5). 



2.2. Data Reduction 

Raw data from the EPIC-PN instrument were repro- 
cessed using Science Analysis Software (SAS) version 6.0 
to produce calibrated event files. Only events with the 
single or double pixels (PATTERN <4) were used in the 
following analysis. Time intervals of background flaring 
were excluded. Source spectra were extracted from a cir- 
cular region centered at the peak of the X-ray source 
counts. Background spectra were taken from an off- 
source region in the same field-of-view and subtracted 
from the source spectra. Appropriate response and an- 
cillary files were created using the RMFGEN and AR- 
FGEN tools in SAS, respectively. All the spectra were 
binned so that each bin contains at least 30 (or more 
for bright sources) counts in order to apply a x 2 mm - 
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TABLE 1 
The Sample 



PG Alternate Name z log (Mbh/Mq) Lbol/-£<Edd Group a Note b Reference 



0003+199 


Mrk 335 


0.026 


7.16 


0.58 


2 




1, 2 


0007+106 


Mrk 1501 


0.089 


8.69 


0.10 


4 


RL 


3, 4 


0050+124 


I Zw 1 


0.061 


7.29 


2.51 


1 




2, 3, 4, 5 


0157+001 


Mrk 1014 


0.163 


8.21 


1.11 


1 




2, 3, 4 


0804+761 




0.100 


8.41 


0.45 


2 




3, 4, 5, 6 


0844+349 


Ton 951 


0.064 


7.80 


0.33 


3 




2, 3, 4, 7, 8, 9 


0947+396 




0.206 


8.67 


0.26 


3 




2, 3, 4, 5, 10 


0953+414 




0.234 


8.69 


0.66 


2 




2, 3, 4, 5, 10 


1004+130 


PKS 1004+13 


0.240 


9.43 


0.09 


4 


RL, BAL 


11 


1048+342 


3C 246 


0.167 


8.36 


0.35 


3 




2, 3, 4, 5, 10 


1100+772 


3C 249.1 


0.312 


9.43 


0.12 


4 


RL 


3, 4 


1114+445 




0.144 


8.57 


0.13 


4 




3, 4, 5, 10, 12 


1115+407 




0.154 


7.61 


1.91 


1 




2, 3, 4, 5, 10 


1116+215 


Ton 1388 


0.177 


8.64 


0.68 


2 




2, 3, 4, 5, 10 


1202+281 


GQ Comae 


0.165 


8.54 


0.18 


3 




2, 3, 4, 5, 10 


1211+143 




0.081 


7.95 


0.79 


2 




2, 3, 4, 13, 14 


1216+069 




0.331 


9.40 


0.18 


3 




3, 4, 10 


1244+026 




0.048 


6.39 


4.57 


1 




2, 3, 4, 5 


1302-102 


PKS 1302-102 


0.27 


9.04 


0.48 


2 


RL 




1307+085 




0.155 


8.87 


0.10 


4 




2, 3, 4, 5 


1309+355 


Ton 1565 


0.184 


8.37 


0.32 


3 


RL 


2, 3, 4, 5, 10, 12 


1322+659 




0.168 


8.22 


0.49 


2 




2, 3, 4, 5, 10 


1352+183 




0.152 


8.37 


0.37 


3 




2, 3, 4, 5, 10 


1402+261 


Ton 182 


0.164 


7.94 


1.78 


1 




2, 3, 4, 5, 10, 15 


1404+226 




0.098 


6.80 


2.63 


1 




2, 3, 4, 16, 17 


1411+442 




0.090 


7.98 


0.28 


3 




3, 4, 10, 18 


1415+451 




0.114 


7.87 


0.38 


3 




3, 4, 10 


1425+267 


Ton 202 


0.366 


9.75 


0.05 


4 


RL 


10, 19 


1426+015 


Mrk 1383 


0.087 


8.98 


0.08 


4 




5, 6 


1427+480 




0.221 


8.10 


0.71 


2 




2, 3, 4, 5, 10 


1440+356 


Mrk 478 


0.079 


7.42 


1.91 


1 




2, 3, 4, 5, 10 


1444+407 




0.267 


8.36 


0.94 


2 




2, 3, 4, 10 


1448+273 




0.065 


6.95 


1.86 


1 






1501+106 


Mrk 841 


0.036 


8.39 


0.07 


4 




2, 3, 4 


1512+370 


4C +37.43 


0.371 


9.50 


0.13 


4 


RL 


5, 10 


1535+547 


Mrk 486 


0.039 


7.02 


0.41 


3 




20 


1543+489 




0.400 


8.13 


3.39 


1 




10 


1613+658 


Mrk 876 


0.129 


9.12 


0.06 


4 




3, 4, 5, 6 


1626+554 




0.133 


8.43 


0.19 


3 




3, 4, 5, 10 


2112+059 




0.466 


9.27 


0.56 


2 


BAL 


20 


2130+099 


Mrk 1513 


0.063 


7.82 


0.65 


2 






2214+139 


Mrk 304 


0.066 


8.38 


0.06 


4 




3, 4, 18 


2233+134 




0.325 


8.15 


2.29 


1 







References. — (1) Gondoin ct al. 2002, (2) Crummy ot al. 2006, (3) Piconcelli ct al. 2005, (4) Jimenez-Bailon et 
al. 2005, (5) Porquet ot al. 2004, (6) Pago ot al. 2004b, (7) Brinkmann ot al. 2006, (8) Brinkmann ot al. 2003, (9) 
Pounds ot al. 2003, (10) Brocksopp ot al. 2006, (11) Miller ct al. 2006, (12) Ashton et al. 2004, (13) Pounds ct al. 
2003, (14) Pounds & Page 2006, (15) Reeves et al. 2004, (16) Crummy ot al. 2005, (17) Dasgupta ot al. 2005, (18) 
Brinkmann et al. 2004, (19) Miniutti & Fabian 2006, (20) Schartel et al. 2005. 

a Group divided by Eddington ratio: (1)1.00 < L bol /L Edd , (2)0.40 < L bo i/£Edd < 1-00, (3)0.17 < L bo i/£Edd < 
0.40, (4)L bol /L Edd < 0.17. 

b RL: Radio-loud quasars (Boroson & Green 1992), BAL: Broad absorption-line quasars (Weymann et al. 1991; Brandt 
et al. 2000). 



imization technique. Spectral fits were performed with 
XSPEC version 12.2. 

3. RESULTS 

3.1. Spectra of Individual Objects 

The spectra of most objects are not well represented by 
a simple power law; soft excess emission and an emission 
line-like feature around 6 keV are also often seen. There- 
fore, we considered a model consisting of a power law, a 
blackbody component to represent the soft excess emis- 
sion below ~ 2 keV, and a Gaussian emission line around 
6 keV, all modified by Galactic absorption. We used 
a two blackbody model for the soft excess if the single 
blackbody did not represent the data well. The hydrogen 
column densities of the Galactic absorption are left free. 



Some objects show absorption features around 0.65-0.95 
keV. We multiplied one or two edges to the above model 
for such cases. 

The spectra of 39 objects are well reproduced by this 
model. We obtained hydrogen column densities of Nn ~ 
(0.003-0.26) xl0 22 cm- 2 , photon indices V w 1.5 to 2.5, 
and temperatures of the blackbody kT « 0.07 — 0.23 keV. 
The column densities are consistent with the Galac- 
tic value derived from the HI map by Dickey & Lock- 
man (1990). Except for a few cases (see Appendix), 
our results for individual objects are in good agreement 
with the independent analysis published in the refer- 
ences given in Table 1. The spectra of the remaining 
four objects (PG 1004+130, 1411+442, 1535+547, and 
2214+139) show broad excess emission at energies above 
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TABLE 2 
Observation Log 



PG 


Obs. Date 


Exposure 
Cks) 


Count rate a 
(count s ^ ) 


Z — 1U KC V 

(10 — ^ergs cm - ^) 


Lo 1 n h-a\r c 

Z — 1U KC V 

(10 44 ergss _1 ) 


Note d 


0003+199 


2000 Dec 25 


28.5 


19.45 


13.75 


0.20 


1 


0007+106 


2000 Jul 03 


10.2 


3.23 


7.01 


1.28 




0050+124 


2002 Jun 22 


18.3 


7.95 


8.17 


0.71 


1 


0157+001 


2000 Jul 29 


4.4 


1.10 


0.91 


0.65 




0804+761 


2000 Nov 04 


0.5 


11.22 


9.94 


2.45 




0844+349 


2000 Nov 04 


10.3 


6.45 


4.87 


0.46 




0947+396 


2001 Nov 03 


17.6 


1.76 


1.78 


2.05 




0953+414 


2001 Nov 22 


10.8 


3.92 


2.94 


4.65 


1 


1004+130 


2003 May 04 


18.2 


0.10 


0.32 


0.48 




1048+342 


2002 May 13 


21.5 


1.17 


1.30 


0.94 




1100+772 


2002 Nov 01 


6.8 


2.02 


3.62 


10.16 




1114+445 


2002 May 14 


35.0 


0.53 


2.28 


1.12 


1 


1115+407 


2002 May 17 


14.4 


2.20 


1.19 


0.75 




1116+215 


2001 Dec 02 


5.5 


5.05 


3.19 


2.77 




1202+281 


2002 May 30 


12.2 


2.64 


3.58 


2.46 


1 


1211+143 


2001 Jun 15 


48.5 


3.33 


3.05 


0.46 


2 


1216+069 


2002 Jul 05 


5.4 


0.82 


1.12 


3.63 




1244+026 


2001 Jun 17 


4.6 


6.15 


2.31 


0.12 




1302-102 


2002 Jan 03 


11.5 


0.61 


3.05 


6.44 




1307+085 


2002 Jun 13 


9.9 


0.80 


1.96 


1.13 


1 


1309+355 


2002 Jun 10 


23.6 


0.44 


0.70 


0.61 


1 


1322+659 


2002 May 11 


8.0 


2.41 


1.27 


0.98 


1 


1352+183 


2002 Jul 20 


8.9 


2.30 


1.87 


1.11 


1 


1402+261 


2002 Jan 27 


9.1 


3.17 


1.80 


1.32 




1404+226 


2001 Jun 18 


9.1 


0.59 


0.13 


0.03 


2 


1411+442 


2002 Jul 10 


21.6 


0.12 


0.46 


0.08 


1 


1415+451 


2002 Dec 08 


20.4 


1.43 


1.06 


0.34 




1425+267 


2002 Jul 28 


29.6 


0.67 


1.68 


6.29 


1 


1426+015 


2000 Jul 28 


2.2 


9.39 


8.48 


1.50 




1427+480 


2002 May 31 


32.5 


1.15 


1.05 


1.44 




1440+356 


2003 Jan 04 


10.8 


5.58 


2.63 


0.39 




1444+407 


2002 Aug 11 


16.7 


0.96 


0.45 


1.05 




1448+273 


2003 Feb 08 


18.0 


4.23 


2.07 


0.20 




1501+106 


2001 Jan 13 


7.6 


18.31 


15.00 


0.44 


1 


1512+370 


2002 Aug 25 


15.4 


1.54 


1.87 


7.98 


1 


1535+547 


2002 Nov 03 


14.4 


0.16 


1.33 


0.04 




1543+489 


2003 Feb 08 


3.7 


0.35 


0.13 


0.87 




1613+658 


2001 Aug 29 


2.4 


4.46 


4.98 


2.05 


1 


1626+554 


2002 May 05 


4.8 


2.84 


2.86 


1.28 




2112+059 


2003 May 14 


7.4 


0.13 


0.34 


2.24 




2130+099 


2003 May 16 


23.3 


2.13 


3.63 


0.32 


1 


2214+139 


2002 May 12 


6.2 


0.55 


3.32 


0.31 


1 


2233+134 


2003 May 28 


6.2 


0.59 


0.49 


1.69 


1 



Count rate in 0.3—15 kcV. 

Flux in 2-10 kcV, not corrected for absorption. 
Luminosity in 2—10 keV, not corrected for absorption. 
Signature of ionized absorber. 1: one edge; 2: two edges. 
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Fig. 1. — Distribution of (a) redshifts, (b) central black hole masses, and (c) Eddington ratios. 
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~ 2-3 keV. We added an absorbed power-law component 
to the model used above to represent this component. 
This model provides a good description of the spectra. 
The hydrogen column densities of the absorbed power 
law are ~ 10 23 cm~ 2 for all the cases. Details of the spec- 
tral fits to these four quasars are given in Appendix A. 
The spectral parameters for the other components are 
similar to those obtained for the 39 objects without an 
additional absorbed power law. 

An absorption feature around 0.65-0.95 keV, which 
can be well fitted by one or two edges, is seen in 20 out 
of the 43 objects. The majority of the sources (18/20) 
require only a single edge, with a mean edge energy of 
E = 0.69 keV and an optical depth of r = 0.56. The 
absorption features in the two remaining objects are rep- 
resented by two edges: for PG 1211+143, the edge en- 
ergies are's = 0.75l°;° 2 /0.93t°;° 2 keV and the optical 

depths are r = 0.36^:^/0.31^;^; for PG 1404+226, 
the prameters are E = 0.74t°;J 2 /0.94l°;^ keV and r = 
0.25lg;2i/0.89lo;37). Similar features arc often seen in 
Seyfert Is and quasars observed with the energy resolu- 
tion of CCD. High-resolution spectroscopy of such fea- 
tures with grating spectrometers onboard XMM-Newton 
and Chandra have shown that they consist of many ab- 
sorption lines of ionized species such as He-like and hy- 
drogenic O, Ne, and Mg. The observed edge energies are 
in agreement with such ionized absorbers observed with 
CCD resolution (e.g., George et al. 1998, 2000). 

The fluxes and luminosities obtained for the best-fit 
model, as well as the detected count rates, are summa- 
rized in Table 2. Neither the fluxes nor the luminosities 
are corrected for absorption. 

3.2. Co-added Spectra 

We derived a co-added rest-frame spectrum of the 43 
quasars in order to examine the average shape of the Fe 
K emission line. Since each object has a different red- 
shift and thus a different observed line energy, the de- 
tector response, which is energy dependent, is different 
for each object. As the co-added rest-frame spectrum is 
a summation of spectra taken with different energy res- 
olution, this situation complicates quantitative analysis. 
Therefore, the spectrum derived here is for presentation 
purpose only, and quantitative analysis by simultaneous 
fitting is done in the next subsection. 

The co-added spectrum was constructed as follows. 
First, ratios of the data to the best-fit continuum model 
were calculated. The ratios were then multiplied by the 
unfolded best-fit continuum model. The resulting spec- 
trum is composed of an unfolded continuum and a folded 
emission- line component. The energy scales of both the 
unfolded continuum spectrum and the spectrum consist- 
ing of the unfolded continuum plus the folded line were 
shifted to the source redshift. We selected bin widths of 
0.25 keV for E < 8 keV and 2.5 keV for E > 8 keV in this 
process, and used a Monte Carlo method to redistribute 
the observed events into the new spectral bins. The 43 
spectra with the folded line are then co-added. The un- 
folded continuum model spectra were also co-added. Fi- 
nally, the ratio of the two spectra is calculated, as shown 
in Figure 2. A prominent narrow core at ~6.4keV and a 
very weak low-energy wing are seen. 



3.3. Simultaneous Fit to All Spectra 

In contrast to previous studies, we employ simultane- 
ous fits rather than fits to stacked spectra. Since the Fe 
K line in each quasar has a different peak energy in the 
observed frame according to their redshifts, a technique 
that properly treats the energy-dependent detector re- 
sponse is necessary. We fitted all spectra simultaneously 
in order to determine mean Fe K line parameters. The 
spectral parameters of the best-fit continuum models for 
each object were used as initial parameters in the spec- 
tral fits. All the continuum parameters were left free, 
while the parameters for the edge model are fixed at the 
best-fit values for individual fits. A Gaussian model with 
common peak energy and width is added to the contin- 
uum model. The normalizations of the Gaussian com- 
ponent were left free individually. The fitting results are 
shown in Table 3. The best- fit peak energy and the dis- 
persion of the Gaussian line are E = 6. 48^0 04 keV and 

a = 0.36^;^ keV, respectively. The mean of the dis- 
tribution of the equivalent widths, shown in Figure 3, is 
248 eV, with a standard deviation of 168 eV. 

A relativistic line model (diskline model in XSPEC) 
instead of a Gaussian was also examined. We fixed sev- 
eral parameters; the peak energy at 6.4 keV, the outer 
radius i? ou t at 500 R Sl where R s is the Schwarzschild ra- 
dius, the emissivity index (3 at —2, and the inclination 
angle of the disk i at 30 degrees. The free parameters 
are the inner radius R{ n and the intensity of the line. 
Again, the parameters in the continuum model were ini- 
tially set to the best-fit values in the individual fits and 
left free except for the edge component. The results are 
summarized in Table 3. The best-fit inner radius i?; n is 
3.0lg Qi? s . The lower boundary is fixed at 3i? s , which is 
the minimum value allowed in the diskline model as- 
suming a Schwarzschild black hole. The mean and stan- 
dard deviation of the equivalent widths are 210 eV and 
133 eV, respectively. 

The line profile in Figure 2 shows a prominent narrow 
core and a very weak low-energy wing. Such a profile 
may be indicative of the presence of multiple line compo- 
nents, and a single Gaussian or a diskline profile may not 
be a good representation. Therefore, we examined a two- 
component line model consisting of either two Gaussians 
or a combination of a Gaussian and a diskline component. 
With regards to the double-Gaussian model, although 
the x 2 improved by 103 for 3 additional parameters, not 
all the line parameters were tightly constrained. Thus, 
we first fixed the continuum component to obtain a rough 
value of the Gaussian widths, and then the widths were 
fixed at these values to constrain the equivalent widths 
of the lines. All the parameters except for the widths 
were left free in the latter fit. Note that the normaliza- 
tions of the lines were determined for each object. The 
resulting peak energies of the broad and narrow lines are 
5.06 keV and 6.52 keV, respectively (see Table 3). The 
mean equivalent width for the narrow component is 261 
eV, and for the broad component it is 150 eV, although 
the standard deviations of the distributions are large (171 
and 181 eV, respectively). The distributions suggest that 
the narrow component is commonly seen and that the 
broad line is present in a fraction of the sample. 

The Gaussian plus diskline model was then examined. 
We fixed the peak energy and the dispersion of the Gaus- 
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sian to the best-fit parameters of the narrow component 
in the double-Gaussian fit. The fixed parameters for the 
diskline component were the same as those in the sin- 
gle cliskline fit presented above. The free parameters are 
the inner radius and the normalization of the diskline, 
and the parameters for the Gaussian and continuum. As 
summarized in Table 3, the best-fit inner radius is i?; n = 
3.0±l% Rs, and the mean equivalent width of the diskline 
is EW=102eV. The large standard deviation (181 cV) 
again suggests that the diskline is required in only a frac- 
tion of the sample. 



Rest Energy (keV) 

Fig. 2. — Data-to-model ratio for co-added spectrum of 43 PG 
quasars. Vertical dotted line is at 6.4 keV. 




Fig. 3. — Distribution of equivalent widths of Fc K line in 43 PG 
quasars. A single Gaussian model is assumed. 



3.4. Eddington Ratio Dependence of the Fe Line Profile 

We now turn to the Eddington ratio dependence of 
the Fe line profile. The 43 quasars were divided into 
four groups according to the Eddington ratio as follows: 
(1)1.00 < L hol /L Edd , (2)0.40 < L hol /L Edd < 1.00, 
(3)0.18 < L bol /L Edd < 0.40, and (4) L hoi /L Edd < 0.18 
(see Table 1) . The spectra in each group were co- added as 
in Section 3.3 for the purposes of presentation; Figure 4 
shows the resulting data-to-continuum model ratios. As 
the Eddington ratio increases, the width of the Fe K line 
appear to becomes broader and its peak energy rises. 

Following the procedure outlined in Section 3.3, we 
fitted the spectra in each group simultaneously to char- 
acterize the Fe K line quantitatively. The following four 



models were examined for the Fe K line component: sin- 
gle Gaussian, single diskline, double Gaussians, and a 
Gaussian plus diskline model. The results of the spec- 
tral fits are shown in Table 4, and the Eddington ratio 
dependence of the peak energy and the width of the line 
obtained from the single Gaussian model fits are shown in 
Figure 5. We again find that the width and peak energy 
of the line become broader and higher as the Edding- 
ton ratio increases. The peak energy gradually increases 
from 6.37 keV in group 4 to 6.77 keV in group 1. The 
width is relatively narrow (a = 0.15 keV) in group 4, but 
definitely resolved (cinst = 0.07 keV), and becomes wider 
(<r = 0.68 keV in group 1) as L bo \/L Edd increases. The 
distributions of the equivalent widths for each group ob- 
tained from the Gaussian fits are shown in Figure 6; the 
mean equivalent widths in each group are 280, 276, 287, 
and 131 eV, respectively. 

A diskline model was examined next. We fixed several 
diskline parameters. The line center energy was chosen 
to be 6.4 keV for groups 3 and 4 and 6.7 keV for groups 
1 and 2, for which the Gaussian fits suggest that the line 
is from ionized Fe. The outer radius i? ut, the emissivity 
index [3, and the inclination angle i were fixed at 500 Rs, 
—2, and 30 degrees, respectively. The free parameters 
are the inner radius i?i n and the normalization. The 
parameters for the continuum model were also left free, 
except for those of the edge component. The results are 
shown in Table 4. The resulting values of \ 2 ar e similar 
to or slightly worse than those for the Gaussian fits. The 
best-fit inner radius is i?; n « 3i?s, which is the last stable 
orbit of a Schwartzschild black hole. 

We also examined a two-component line model since a 
strong narrow core and a broad wing are visible in groups 
3 and 4, and possibly in group 2. A double-Gaussian 
model was first tested in the same way as in Section 3.3; 
the results are shown in Table 4. The resulting x 2 values 
are similar to or only slightly better than those of the 
single- Gaussian fit, indicating that the additional broad 
line is statistically not required. The distributions of the 
EW are broad and indicate that some objects may have 
a broad component. 

A Gaussian plus diskline model was then used in or- 
der to constrain the equivalent width of the diskline, in 
which the inner radius is fixed at 3Rs- The peak energy 
and the dispersion of the Gaussian were fixed at the best- 
fit parameters for the narrow component in the double- 
Gaussian fit. The inner radius of the disk was also fixed 
at 3i?s, and the other parameters in the diskline model 
were fixed to the same values as in the single-diskline fits 
described above. Only the normalizations of the Gaus- 
sian and diskline were left free. The parameters for the 
continuum model were treated in a similar way as in 
Section 3.3. The results are summarized in Table 4. The 
resulting values of x 2 are similar to or only slightly better 
than those for the single-Gaussian or diskline fits. The 
mean equivalent widths of the diskline in each group are 
99, 176, 132, and 67 eV, respectively, although the stan- 
dard deviations as large. The large standard deviations 
again suggest that the broad diskline is not ubiquitous. 

4. DISCUSSION 

4.1. Origin of the Eddington Ratio Dependence 

We examined mean Fe K line shapes for the four groups 
sorted by the Eddington ratios, and found that the peak 
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TABLE 3 
Results of Spectral Fit. 



Gaussian Disklinc 



Model 


Energy 
(kcV) 


(keV) 


EW a 
(eV) 


Energy 
(kcV) 


Rin 

(R B ) 


Rout 

(Rs) 


P 


(deg) 


EW a 
(eV) 


X 2 /d.o.f. 


Gaussian 


6-48l u :S 


u - oo -0.08 


248(168) 














8958/8129 


Diskline 








6.4f 




500f 


-2f 


30f 


210(133) 


8975/8130 


Gaussian+Gaussian 


5.06lg;g 


0.61f 


150(181) 


















°- oz -0.04 


0.34f 


261(171) 














8853/8086 


Gaussian+Diskline 


6.52f 


0.34f 


155(181) 


6.4f 


, n+2.9 
^■"-O.Op 


500f 


-2f 


30f 


102(122) 


8912/8087 



a Numbcrs in parentheses are standard deviation of distribution of 
best-fit values. 

"p" denotes pegged parameter. 
c "f " denotes fixed parameter. 



TABLE 4 

Results of Spectral Fit in Each Eddington Ratio Group. 



Gaussian Disklinc 



Group 


Model 


Energy 
(keV) 


a 

(keV) 


EW a 
(eV) 


Energy 
(keV) 


Rin 

(Rs) 


Rout 

(Rs) 


P 


(deg) 


EW a 
(eV) 


*7d.o.f. 


1 


Gaussian 


D - ' ' -0.30 


0-68^0.24 
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2308/2237 




Diskline 








6.7f 


r, n +20.o 

aM -0.0p 


500f 


-2f 


30f 


202(176) 


2316/2238 




Gaussian+Gaussian 


6 58+ 014 


0.42f 


119(152) 


















Gaussian+Diskline 


- 4 ~+0.58 

-1.44 
7.47f 


0.48f 
0.48f 
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208(281) 


6.7f 


3.0f 


500f 


-2f 


30f 


99(158) 


2305/2248 
2304/2250 


2 


Gaussian 


6 54+ U uy 

D -°^-0.08 


U - JO -0.07 


276(160) 














2636/2197 




Diskline 








6.7f 


3 0+ 19 


500f 


-2f 


30f 


285(171) 


2662/2198 




Gaussian+Gaussian 


5 02+ 30 


0.38f 


90(137) 


















Gaussian+Diskline 


6.52f 


0.24f 
0.24f 


183(96) 
68(80) 


6.7f 


3.0f 


500f 


-2f 


30f 


176(203) 


2634/2218 
2646/2220 


3 


Gaussian 


6 45+ u (Jy 


U -" J ' -0.12 


287(179) 














1796/1774 




Diskline 








6.4f 


o o+S.O 


500f 


-2f 


30f 


259(149) 


1801/1775 




Gaussian+Gaussian 


c; t:c;+0.46 
°-°°-0.S7 


0.48f 


96(123) 


















Gaussian+Diskline 


fi AO+ - 09 
b - 4y -0.08 

6.49f 


0.25f 
0.25f 


173(124) 
83(114) 


6.4f 


3.0f 


500f 


-2f 


30f 


132 (132) 


1800/1798 
1803/1800 


4 


Gaussian 


6 37 +U.U& 

— .04 


u - ±J -0.02 


131(59) 














2184/1935 




Diskline 








6.4f 


o n +24.2 


500f 


-2f 


30f 


179(68) 


2186/1936 




Gaussian+Gaussian 


5 90+ - 37 


0.70f 


104(108) 


















Gaussian+Diskline 


fi o 7 +0.04 
— .03 

6.37f 


0.09f 
0.09f 


82 (70) 
75(76) 


6.4f 


3.0f 


500f 


-2f 


30f 


67(73) 


2173/1957 
2180/1959 



a Numbcrs in parentheses are standard deviation of distribution of 
best-fit values. 

"p" denotes pegged parameter. 
c "f " denotes fixed parameter. 

energy, width, and equivalent width become higher (6.4 
to 6.8keV), broader (0.1 to 0.7keV), and larger (130 to 
280 eV), respectively, as the Eddington ratio increases. 
In this section, we discuss possible reasons of the ob- 
served Eddington ratio dependence of the Fc line profile. 

The line center energy of the Fe K line is governed by a 
combination of the Dopplcr effect caused by the orbital 
motion, gravitational redshift, and the ionization state 
of the emitter. An Fe line from a rotating disk around 
a black hole shows a double-peaked profile with a peak 
energy that depends on the inclination angle of the disk. 
Since each of the four groups we used contains from 10 
to 11 objects, it is unlikely that any one of the groups 
has a mean inclination angle considerably different from 
those of other groups. Therefore, the different line center 
energies obtained for the four groups probably reflect dif- 
ferent ionization states rather than different mean disk 
inclination angles. The Eddington ratio dependence of 



the peak energy thus strongly suggests that the ioniza- 
tion state of the accretion disk depends on the Edding- 
ton ratio. Many aspects of an ionized disk have been 
studied theoretically over the last decade. Calculations 
have been performed for constant-density atmospheres 
(e.g., Matt et al. 1993; Ross & Fabian 1993; Zycki ct 
al. 1994), as well as for atmospheres in hydrostatic equi- 
librium (Nayakshin et al. 2000; Ballantyne et al. 2001). 
These authors calculated reflection spectra from an ion- 
ized disk for different conditions. In the lowest ionization 
regime, the most prominent feature is the cold iron line at 
6.4 keV. As the mass accretion rate goes up, the ionized 
stage becomes higher, and the Fe K line is dominated by 
Fe XXV and XXVI, and K lines from the lighter elements 
emerge in the 0.3 — 3keV band. At very high accretion 
rates, the surface of the disk is highly ionized so that the 
only noticeable line is a Compton-broadened Fe K line 
peaking at ^7keV. The Eddington ratio dependence of 
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Fig. 4. — Eddington ratio dependence of the Fe line profile for 
each of the four groups. The vertical dotted line is at 6.4 keV. 



the peak energy and the width of the Fe K line we ob- 
served in the PG quasar sample is in good agreement 
with the predictions of ionized disk models. 

The large EW in the high-Eddington ratio groups is 
also likely to be related to the ionization state of the ac- 
cretion disk. Calculations of the reflected X-rays from an 
ionized disk have shown that the EW of the Fe K line can 
be higher than in the case of a disk in a low-ionization 
state. The resonant trapping effect plays a role in the 
intermediate-ionization state higher than Fe XVIII. The 
EW increases again from Fe XX, reaching a maximum 
value of 250 - 300 eV (Matt et al. 1993; Ross & Fabian 
1993; Zycki & Czerny 1994). The observed line center 
energies and EWs for the three groups with a largest Ed- 
dington ratios are consistent with the latter expectation, 




-0.5 o 

Fig. 5. — Eddington ratio dependence of the peak energy and 
the width of the line. The horizontal dotted line on the top panel 
is at 6.4 keV, and that on the bottom panel is at 0.07 keV, the 
instrumental resolution of XMM-Newton. 




200 300 
EW (eV) 

Fig. 6. — Distribution of Fe line equivalent width in each Ed- 
dington ratio group. 



although the scatter of the actual EW distributions is 
large. The reduction of EW caused by resonant trapping 
is not clearly seen in our data. This may be due to the 
coarse binning of the Eddington ratio we have adopted. 
Future analysis of a larger sample may be able to see the 
ionized disk in an intermediate-ionization state. 

The Eddington ratio dependence of the line width also 
suggests a change of the geometrical structure of the disk, 
if the velocity width indeed reflects emission from the in- 
ner part of an accretion disk broadened by Doppler ef- 
fects and gravitational redshift. Although the line widths 
obtained from the single-Gaussian fits are consistent with 
systematic variations in the inner disk radius i?i n , such a 
trend is not clear in the fits obtained using the disklinc 
model. It is possible that the broader line widths at high 
Eddington ratios may in part be accounted for by a blend 
of lines from Fe in different ionization states. 

4.2. Comparison with Average Fe Line Profile of AGNs 
in Lockman Hole 

Strcblyanska ct al. (2005) derived an average Fe line 
profile of 53 type 1 AGN detected in a 770 ks observation 
of the Lockman Hole with XMM-Newton. Their average 
spectrum shows a broad line-like feature (a = 0.69 keV) 
peaking at 6.4 keV with an equivalent width of EW = 
420 eV. The profile is asymmetric, with a significant wing 
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seen toward low energies. The line width (a — 0.36 keV) 
we obtained for the PG quasar sample is significantly nar- 
rower, and the equivalent width (EW«250eV) smaller, 
than that of Streblyanska et al.'s profile. We discuss pos- 
sible causes for the difference. 

4.2.1. Contribution of Absorbed Quasars 

Appropriate modeling of the continuum is essential to 
derive a correct Fe line profile. A single power-law model 
was assumed in the analysis of the AGNs in the Lockman 
Hole. Since the spectra of many AGNs show a signature 
of absorption by ionized and/or cold matter, a simple 
power law may not be a good approximation of the un- 
derlying continuum. Indeed, four quasars in our sam- 
ple are absorbed by a large column density (N-r w 10 23 
cm -2 ). If there are some highly absorbed objects in the 
Lockman Hole sample, using a single power-law model 
for the continuum may mimic a broad-line profile. 

The flux of a highly absorbed power law with a column 
density of 10 23 cm~ 2 at 5 keV is about 50% of that of an 
unabsorbed power law. If the fraction of highly absorbed 
quasars is 1/10, their contribution to the low-energy side 
of the Fe K line profile around 5 keV is less than 5%. The 
contribution of the absorbed continuum depends on the 
fraction of absorbed quasars. Mateos et al. (2005) mea- 
sured the absorbed fraction for their sample of Lockman 
Hole AGNs with sufficient X-ray counts. They found 
that seven out of 46 type 1 AGNs show a highly ab- 
sorbed spectrum. This fraction is similar to that in our 
PG quasar sample. Note, however, that Streblyanska et 
al.'s sample may contain more highly absorbed objects 
because their flux limit is lower than that used by Ma- 
teos et al. (2005) and X-ray surveys have shown that 
objects with lower flux tend to have a harder spectrum. 
If the absorbed fraction is 10% or higher, rather than 
5%, in the Lockman hole sample, a considerable fraction 
of the red wing might be explained by the contribution 
of highly absorbed emission. 

4.2.2. Distribution of Luminosities 

The physical state of an accretion disk may depend on 
luminosity. If the mean luminosity of our sample is fairly 
different from that of the Lockman Hole AGNs, it may 
introduce a difference in the Fe K line profile. In order to 
examine this possibility, we compared the distributions 
of the 2 — lOkeV luminosity between the two samples 
(Fig. 7). The luminosities of the Lockman Hole AGNs 
were taken from Mainieri et al. (2002). They analyzed 
39 type 1 AGNs in the XMM-Newton observation of the 
Lockman Hole field. Most of the 53 type 1 AGNs in Stre- 
blyanska et al. (2005) have been analyzed by Mainieri 
et al. (2002). The two distributions are very similar, 
and their mean luminosities are essentially identical: the 
mean luminosity of our sample is 10 43 ' 98 ergss -1 , while 
that of the Lockman Hole AGNs is 10 44 00 ergs s" 1 . We 
conclude that luminosity differences is not the main cause 
of the difference of the line profiles. 

4.2.3. Artificial Distortion of Spectra 

Yaqoob (2005) pointed out a spurious effect in adding 
spectra of faint sources. He showed that conventional 
methods for averaging low signal-to-noise ratio X-ray 
spectra result in a weakening of emission lines, a broad 



- PG quasars 



Lockman Hole 
type 1 AGNs 



43 43 44 45 46 

lo « L 2-10 t.v 

Fig. 7. — Distribution of luminosities for the PG quasars in our 
sample and for the AGNs in the Lockman Hole. Mean luminosities 
are Lpg=10 43 ' 98 ergss -1 and Llh= 10 44 00 ergss -1 , respectively. 



dip in the continuum above an emission line, and a spec- 
tral hardening at the highest energies, if sources at dif- 
ferent redshifts are stacked. Yaqoob (2005) simulated 
200 spectra from sources with redshifts 0.5-2.5 and 2- 
10 keV fluxes 6-60xl0~ 14 ergs cm~ 2 s -1 and found an 
artifact wing at energies below the peak of the Fe K 
line. This feature occurs at the level of about 5%-10% 
at 4-6 keV. This kind of spectral distortion may affect 
the Fe line shape in the analysis of Streblyanska et al. 
(2005) since their sample has a relatively low average flux 
(^2-iokoV ~ 10~ 140 ergs s" 1 cm -2 ), compared to our 
PG quasar sample (i*2-iokcV ~ 10~ 115 ergs s _1 cm -2 ). 

5. CONCLUSIONS 

We presented the average Fe K line profile of 43 PG 
quasars observed by XMM-Newton. Our findings are 
summarized as follows. 

• The spectra of 39 objects are well reproduced by a 
canonical model consisting of a power law, a black- 
body, an edge in the case of objects with evidence 
for warm absorption, and a Gaussian, all modified 
by Galactic absorption, while that of the remaining 
four objects require additional absorption with a 
column density of 10 23 cm~ 2 to the power-law com- 
ponent. 

• An average spectrum of the 43 quasars shows a 
prominent narrow Fe K line (E — 6.48 keV, a — 
0.36 keV, and EW = 248 cV). 

• The average Fe K line is well represented by dou- 
ble Gaussians or a combination of a Gaussian and a 
diskline model. The best-fit inner radius of the disk 
is 3i? s (< 5.9i? s ) and the mean EW is -100 eV. 
The large standard deviation of the EW distribu- 
tion indicates that a broad disklike line is not ubiq- 
uitous. 

• We found an Eddington ratio dependence of the Fe 
K line profile. As the Eddington ratio increases, 
the peak energy becomes higher (6.4 to 6.8 keV), 
the width becomes broader (0.1 to 0.7keV), and 
the EW is larger (130 to 280 eV). The higher en- 
ergy peak and larger EW at high Eddington ratio 
can be explained by the ionization of the accretion 
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TABLE 5 

Spectral Parameters for Highly Absorbed Quasars. 
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kT 


(A r H)z=o a N H h 


Edge/r 




a 


EW 


X 2 /d.o.f. 
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-0.02 
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a Foreground absorption within our Galaxy 
b Absorption within quasar 

disk. The systematic broadening of the line with in- 
creasing Eddington ratio may indicate a systematic 
decrease of the inner radius of the disk, or perhaps 
blending from a combination of Fe lines in different 
ionization states. 

• The line width we obtained is narrower and the 
equivalent width smaller than those derived from 



the stacked spectrum of AGNs in the Lockman 
Hole (a = 0.69 keV, EW = 420 eV). We attribute 
these differences to either an increased contribution 
of absorbed sources in the Lockman Hole compos- 
ite spectrum or to artifacts introduced in averaging 
faint sources. 



APPENDIX 
NOTES ON INDIVIDUAL OBJECTS 
Highly Absorbed Quasars 

We found four highly absorbed quasars in the 43 PG quasars. We show their spectra in Figure 8. They are not well 
reproduced by the canonical model, and an additional highly absorbed component is required. Therefore, we added an 
absorbed power-law component. The fitting results are summarized in Table 5. The average column densities of the 
additional absorption are ~ 10 23 cm~ 2 . Spectra of PG 1411+442 and PG 2214+139 are presented in Brinkmann et 
al. (2004). The results of PG 1411+442 are consistent with ours. The reduced \ 2 value is not good for PG 2214+139 
(xt ~ 1-38) because of non-systemtaic deviations of the data from the model seen over the energy range of 0.5-4 
keV. Miller et al. (2006) and Schartel et al. (2005) found absorbed emission in PG 1004+130 and PG 1535+547, 
respectively. The results obtained by Miller et al. (2006) are also consistent with ours. Schartel et al. (2005) found the 
peak energy of Fe K line at 6.00lg;2i keV and the width of 0.261q; 2 ^ keV, while the line center energy is lower than 5 keV 
in our analysis. Therefore, we fixed it at 6.4 keV and repeated the fits. The best-fit line width is a — 0.00 (< 0.24) keV 
(see Table 5). 

A population of soft X-ray-weak (SXW) quasars, where the soft X-ray flux is ~ 10 — 30 times smaller than in typical 
quasars (Yuan et al. 1998), has been known. This class comprises ~10% of optically selected quasars. Brandt et al. 
(2000) found 10 SXW quasars in the Boroson & Green (1992) sample of 87 PG quasars. All of the highly absorbed 
quasars presented here are consistent with the SXW classification of Brandt et al. (2000). The large absorption seems 
to be the cause of the X-ray weakness at least in these four objects. 

Quasars Newly Analyzed in This Paper 

We summarize the results of our spectral analysis of quasars newly presented in this paper. Their spectra are shown 
in Figure 9. We fitted the spectra with the canonical model described in Section 3.1; results are listed in Table 6. They 
are well reproduced by the model except for one case (PG 1448+273). The soft excess component in PG 1448+273 is 
not explained by a single blackbody component. We added an additional blackbody component and obtained a good 
fit. 

Quasars Reanalyzed in This Paper 

This section shows results of quasars reanalyzed in this paper, for which the best-fit parameters are not consistent 
with previously published results. Our results are summarized in Table 7. Some notes on individual objects are given 
below. 

PG 0003+199.— Gondoin et al. (2002) give the line center of Fe K line at 6.0+0. 14 keV and the width of 1.2+ 
0.2 kcV. We obtained a higher center energy 6.60^;^ keV and a narrower width 0.351q'qq keV. 

PG 0804 + 761. - Jimcncz-Bailon et al. (2005), Porquet et al. (2004), and Page et al. (2004) show the best-fit 
parameters of the line center at 6.38 ± 0.06, 6.671q'37, and 6.62 ± 0.14keV, respectively. The peak energy measured 
by Porquet et al. (2004) and Page et al. (2004b) are consistent with our result (6.57 ± 0.08 keV). 

PG 12U+026.— The best-fit line center energy is 6.13^Q^ 1 2 keV in our analysis. Jimcnez-Bailon et al. (2005) and 
Porquet et al. (2004) obtained the line center energy at 6.65lg;°g keV and 6.66±™? keV, respecrively. 




TABLE 6 

Spectral Parameters for Quasars Newly Analyzed in This Paper 



PG 


r 


kT 


(A r H)z=O a 


Edge/r 




cr 


EW x 2 /d.o.f. 






(keV) 


(10 22 cm- 2 ) 


(keV) 


(keV) 


(keV) 


(eV) 



1302-102 
1448+273 

2130+099 
2233+134 



1.66 
2.25 



+0.10 
0.11 
+0.11 
0.11 



1.65] 



-+0.03 
-0.02 

2 17+ - 19 

z - i ' -0.17 



0.20 
0.12 
0.23 
0.11 
0.19 



-0.02 
-0.04 
+0.01 
0.02 
+0.06 
0.14 
0.00 
-0.00 
-0.03 



0.00 
0.05 



+0.06 
0.00 

+0.01 
0.01 



0.05 
0.0 



+0.01 
-0.01 
k+0.03 



0.66 



+0.02 /0 27 +0.07 



72+ - 02 in 2Q+ - 06 
u - ' z -0.0l/ u - ZM -0.07 



0.89 i 



70.63^ 



6.41 



+0.54 

0.24 



°- 4o -0.12 



O.lf 



0.02 



+0.22 
0.02 



93 



139 



-142 
-93 



+495 
139 



64/63 
260/290 

182/139 
117/96 



Note. — PG 1448+273 requires two blackbody components to represent soft excess. 
a Foreground absorption within our Galaxy 



TABLE 7 

Spectral Parameters for Quasars Reanalyzed in This Paper. 



PG 


r 


kT 
(keV) 


(A r H)z=0 a 
(10 22 cm- 2 ) 


Edgc/r 
(keV) 


(keV) 


cr 

(keV) 


EW 
(eV) 


X 2 /d-o.f. 


0003+199 
0804+761 
1244+026 
1613+658 


o 9C+O.OI 
z - zo -0.01 
9 o n +0.09 
z.ou_ 09 

2 54+ ' 01 

^•y^-o.oi 

, Q9+O.IO 
^"-O.IO 


n ,,+0.00 

15+ - 01 

1 S+ 002 
u - lo -0.03 


03+ tM)() 

02+ ' 01 
u ' uz -0.01 

0-0i±g;gi 

0.04l° :u 2 


0-65±8;gi /0.l6l ;° 2 

n 7 o+0.16 /n 99+O.I9 
u - '°-0.07/ u - zz -0.20 


a cn+0-13 

6.60I .o 2 

D -°' -0.08 
6 12 


S5+ 016 
O.lf 
O.lf 


192+if 
1406l 2 ^ 3 
6421^1 


626/483 
141/138 
197/187 
109/97 



a Foreground absorption within our Galaxy 

PG 1613+658. — A line feature is not clearly seen at ~6.4keV in our spectrum. No detection is consistent with the 
result of Page et al. (2004a) (EW<85 eV), while Page et al. (2004b) reported EW=96 ± 59 eV. 
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